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^ ' Binaries have played a crucial role many times in the history of modern astronomy and 

\l , are doing so again in the rapidly evolving exploration of the Kuiper Belt. The large fraction of 

' transneptunian objects that are binary or multiple, 48 such systems are now known, has been 

an unanticipated windfall. Separations and relative magnitudes measured in discovery images 
give important information on the statistical properties of the binary population that can be 
related to competing models of binary formation. Orbits, derived for 13 systems, provide a 
\ determination of the system mass. Masses can be used to derive densities and albedos when 

■~>^ ■ an independent size measurement is available. Angular momenta and relative sizes of the 

' majority of binaries are consistent with formation by dynamical capture. The small satellites 

of the largest transneptunian objects, in contrast, are more likely formed from collisions. 
^ ■ Correlations of the fraction of binaries with different dynamical populations or with other 

5^ ' physical variables have the potential to constrain models of the origin and evolution of the 

, transneptunian population as a whole. Other means of studying binaries have only begun to be 

exploited, including lightcurve, color, and spectral data. Because of the several channels for 
obtaining unique physical information, it is already clear that binaries will emerge as one of the 
most useful tools for unraveling the many complexities of transneptunian space. 
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HISTORY AND DISCOVERY populations can be related to formation and subsequent evo- 

lutionary and environmental conditions. 



Ever since Herschel noticed, two hundred years ago, that 
gravitationally bound stellar binaries exist, the search for 



Searches for bound systems among the small body pop- 

, . . , ^ „ , , , , , ^ , ,. ulations in the solar system has a long and mostly fruitless 

binaries has followed close on the heels or the discovery .i..,.!. j- i 

■' history that has been summarized in several recent reviews 

of each new class of astronomical obiect. The reasons for ,. ^ , ^^vrvi a? 77 ^nnA ^rvrv^ n- 1 j j 

■' (Merline et al, 2002; Noll, 2004, 2006; Richardson and 

such searches are, of course, eminently practical. Binary or- ,1,7, -,rvrvc\ r. . j- r 

■' ^ ■' Walsh, 2005). But the discovery of the second transnep- 

bits provide determinations of system mass, a fundamental ^ . /t,x,„\ inno i\r\t7 /t/ -77 . 7 ^rvrvi\ 

, . , .... ■ ,-rn , ■ , tunian binary (TNB), 1998 WW31 (Veillet et al., 2001) 
physical quantity that is otherwise difficult or impossible to 



marked the start of a landslide of discovery that shows no 
signs of abating. 



obtain. The utilization of binaries in stellar astronomy has 
enabled countless applications including Eddington's land- 
mark mass-luminosity relation. Likewise, in planetary sci- , , t^. j r^. j. • j.- e r^x. 

; , • , , , , 1-1 Discovery and Characterization of Charon 

ence, bound systems have been extensively exploited; they 

have been used, for example, to determine the masses of r^, 11 ^ ^ r , u n t-xti-> 

, „ , . . „ , The nrst example of what we would now call a TNB 

planets and to make Roemer s nrst determination of the ,. ■ ■ ■ -.-rr , , ■ , , 

^ , ^ ,. , rr,, . . . . . , was discovered during a very different technological epoch 

speed of light. The statistics of binaries in astronomical . • . j j . 1 

^ than the present, prior to the widespread astronomical use 
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Fig. 1. — Images of 2001 FLiss (top-left), (42355) Typhon/Echidna (top-right), 1999 OJ4 (bottom-left), and 2002 GZ31 (bottom-right). 
The images shown are each combinations of four separate 300 s exposures taken with the High Resolution Camera on HST. The dithered 
exposures have been combined using multidrizzle. The images are shown with a linear grayscale normalized to the peak pixel. The 
pixels in the drizzled images are 25 milliarcsec on a side in a non-distorted coordinate frame. Each of the four panels is 1 arcsec square. 
Images are oriented in detector coordinates. 



of CCD arrays. The discovery of Charon (Christy and Har- 
rington, 1978) on photographic plates taken for Pluto as- 
trometry heralded a spectacular flourishing of Pluto science, 
and offered a glimpse of the tremendous potential of TNBs 
to contribute to Kuiper belt science in general. 

Charon's orbit (Christy and Harrington, 1978, 1980) re- 
vealed the system mass, which up until then had been es- 
timated via other methods, with wildly divergent results. 
About the same time, spectroscopy revealed the presence 
of methane on Pluto (Cruikshank et al, 1976) indicating a 
high albedo, small size, and the possible existence of an at- 
mosphere. Observations of occultations of stars by Pluto 
confirmed its small size (Millis et al, 1993) and enabled 
direct detection of the atmosphere (Hubbard et al, 1988). 

The orbit plane of Charon, as viewed from the Earth, 
was oriented edge-on within a few years of Charon's dis- 
covery. This geometry happens only twice during Pluto's 
248 year orbit, so its occurrence just after Charon's dis- 
covery was fortuitous. Mutual events, when Charon (or its 
shadow) passed across the face of Pluto, or Pluto (or its 
shadow) masked the view of Charon, were observable from 
1985 through 1990 (Binzel and Hubbard, 1997). From the 
timing of these events and from the changes in observable 
flux during them, much tighter constraints on the sizes and 
albedos of Pluto and Charon were derived (e.g. Young and 
Binzel, 1994). Mutual events also made it possible to dis- 
tinguish the surface compositions of Pluto and Charon (e.g. 
Buie et ah, 1987; Fink and DiSanti, 1988), by comparing 
reflectance spectra of the two objects blended together with 
spectra of Pluto alone, when Charon was completely hidden 



from view. From subtle variations in flux as Charon blocked 
different regions of Pluto's surface (and vice versa), maps 
of albedo patterns on the faces of the two objects were con- 
structed (e.g. Buieetal., 1997; Young et al, 1999). The mu- 
tual events are now over and will not be repeated during our 
lifetimes, but telescope and detector technology continue 
their advance. For relatively well-separated, bright TNBs 
like Pluto and Charon, it is now feasible to study them as 
separate worlds even without the aid of mutual events (e.g. 
Brown and Calvin, 2000; Buie and Grundy, 2000). Just as 
the discovery of Charon propelled Pluto science forward, 
the recent study of two additional moons of Pluto (Weaver 
et al, 2006) can be expected to give another boost to Pluto 
science, by enabling detailed studies of the dynamics of the 
system (Buie et al., 2006; Lee and Peale, 2006), and pro- 
viding new constraints on formation scenarios (e.g. Canup, 
2005; Ward and Canup, 2006). 

1.2 Discovery of Binaries 

The serendipitous discovery of the second TNB, 1998 
WW31 (Veillet et al, 2001) marked a breakthrough for bi- 
naries in the Kuiper Belt. It immediately provided a context 
for Pluto/Charon as a member of a group of similar systems 
rather than as a unique oddity. The relatively large sep- 
aration and size of the secondary dispelled the notion that 
Kuiper Belt satellites would all be small, faint, and difficult- 
to-resolve collision fragments. Some, at least, were de- 
tectable from the ground with moderately good observing 
conditions as had been foreseen by Toth (1999). The next 
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two discovered TNBs were just such systems (Elliot et ai, 
2001; Kavelaars et ai, 2001). 

The first conscious search for satelhtes of transneptunian 
objects (TNOs) was carried out by M. Brown and C. Trujillo 
using the Space Telescope Imaging Spectrograph (STIS) on 
the Hubble Space Telescope (HST) starting in August 2000 
(Trujillo and Brown, 2002; Brown and Trujillo, 2002). A 
series of large surveys with HST followed producing the 
discovery of most of the known TNBs (e.g. Noll et ai, 
2002a,b,c, 2003; Stephens et ai, 2004; Stephens and Noll, 
2006; Noll et ai, 2006a,b,c,d, e,f; Figure 1). 

Large ground-based surveys - Keck (Schaller and 
Brown, 2003), Deep Ecliptic Survey (DES) followup at 
Magellan (Millis et ai, 2002; Ellliot et ai, 2005 ), and 
the Canada-France Ecliptic Plane Survey (CFEPS; Allen 
et ai, 2006) - have produced a few detections. Though 
significantly less productive because of the limited angular 
resolution possible from the ground compared to HST, the 
sheer number of objects observed by these surveys makes 
them a valuable statistical resource (Kern and Elliot, 2006). 
Both space- and ground-based discoveries are described in 
detail in §2.2. 

Binaries may also be "discovered" theoretically. Agnor 
and Hamilton (2006) have shown that the most likely expla- 
nation for the origin of Neptune's retrograde satellite Triton 
is the capture of one component of a binary that encoun- 
tered the giant planet. The viability of this model is en- 
abled by the paradigm-shifting realization that binaries in 
the transneptunian region are common. 

2. INVENTORY 

Much can be learned about binaries and the environment 
in which they formed from simple accounting. The fraction 
of binaries in the transneptunian population is far higher 
than anyone guessed a decade ago when none were yet rec- 
ognized (except for Pluto/Charon), and considerably higher 
than was thought even four years ago as the first spate of 
discoveries was being made. As the number of binaries 
has climbed, it has become apparent that stating the frac- 
tion of binaries is not a simple task. The fraction in a given 
sample is strongly dependent on a number of observational 
factors, chiefly angular resolution and sensitivity. To add to 
the complexity, TNOs can be divided into dynamical groups 
with possibly differing binary fractions. Thus, the fraction 
of binaries in a particular sample also depends on the mix 
of dynamical classes in the sample. Perversely, perhaps, the 
brightest TNOs, and thus the first to be sampled, belong to 
dynamical classes with a lower overall fraction of sizable 
binary companions. Other dependencies may also help de- 
termine the fraction and nature of binaries and multiples; 
for example, very small, possibly collision-produced com- 
panions appear to be most likely around the largest of the 
TNOs. The current inventory of TNBs, while impressive 
compared to just a few years ago, remains inadequate to 
address all of the questions one would like to ask. 



2.1 Current Inventory of Transneptunian Binaries 

As of February 2007, more than 40 TNO and Centaur bi- 
naries had been announced through the International Astro- 
nomical Union Circulars (lAUC) and/or other publications. 
Additional binaries are not yet documented in a publication. 
All of the binaries of which we are aware are compiled in 
Table 1 . References listed are generally the discovery an- 
nouncement, when available. Values for separation and rel- 
ative magnitude were recalculated for many of the objects 
and supersede earlier published values. The osculating he- 
liocentric orbital parameters a©, e©, and zq are listed as 
well as the dynamical class. For the latter we have followed 
the DES convention (Elliot et ai, 2005), with the resonant 
objects identified by their specific resonance as n:m, where 
n refers to the mean motion of Neptune. Classical objects 
on orbits of low inclination and eccentricity are designated 
C. Classical objects with an integrated average inclination 
i > 5° relative to the invariable plane are denoted by H. 
Objects in the scattered disk are labelled S or X depend- 
ing on whether their Tisserand parameter, T, is less than or 
greater than 3. The Tisserand parameter relative to Neptune 
is defined as T/v = aM/a + 2[(1 — e'^)a/aN]^/^ cos{i) 
where a,e, and i are the heliocentric semimajor axis, eccen- 
tricity, and inclination of the TNO and is the semimajor 
axis of Neptune. Objects in the extended scattered disk, X, 
have the additional requirement of a time-averaged eccen- 
tricity greater than 0.2. Centaurs and Centaur-like0 objects, 
labelled 0, are on unstable, non-resonant, planet-crossing 
orbits and are, therefore, dynamically young. In Table 1 we 
list the objects in three broad dynamical groupings. Classi- 
cal, Scattered, and Resonant. The Classical grouping in- 
cludes all Classical objects regardless of inclination, i.e. 
both Hot and Cold Classicals. The Resonant group includes 
all objects verified to be in mean motion resonances by nu- 
merical integration. The Scattered group includes both Near 
and Extended Scattered objects and the Centaurs. Within 
each group we have ordered the objects by absolute magni- 
tude, Hv- 

In addition to the dynamical class and osculating he- 
liocentric orbital elements of the binaries, we list in Ta- 
ble 1 three additional measurements available for all of the 
known binaries: 1 .) The reported separation at discovery 
(in arcsec) is shown with the error in the final significant 
digit in parentheses. Separations reported without an error 
estimate are shown in italics. As we discuss in more detail 
below in §3.3, the separation at discovery is not an intrinsic 
property of binary orbits, but can be useful for estimating 
the distribution of binary semimajor axes. 2.) The magni- 
tude difference, A,nag, can be used to derive the size ratio 
of the components (with the customary assumption of equal 
albedos). Once again, errors in the final digit are shown in 



^Centaur-like objects are in unstable, non-resonant, giant-planet-crossing 
orbits just like the Centaurs, but have a semimajor axis greater than 30. 1 
AU. There is currently disagreement on what this class of objects should 
be called. Because of their similarly unstable orbits, we follow the DES 
convention and group them with Centaurs. 
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parentheses, and estimated quantities are in italics. 3.) The 
absolute magnitude, Hv, is taken from the Minor Planet 
Center (MPC) and applies to the combined light of the un- 
resolved binary. Better measurements of Hy are available 
for some objects (Romanishin and Tegler, 2005), but for the 
sake of uniformity we use the MPC values for all objects. 
Hv can provide a determination of the size if the albedo is 
known or can be estimated. However, the range of albedo in 
the transneptunian population is large {Grundy et ai, 2005) 
as is the phase behavior {Rabinowitz et al, 2007) making 
any such estimate risky (see also the chapters by Stansberry 
et al. and Belskaya et al). 

2.2 Large Surveys, Observational Limits, and Bias 

Several large surveys have produced the discovery of 
large fractions of the known TNBs. Observational limits 
are, to first order, a function of the telescope and instrument 
used for the observations. This is more easily characterized 
for space-based instruments than for ground-based surveys, 
but approximate limits for the latter can be estimated. 

The largest semi-uniform ground-based survey of TNOs 
that has been systematically searched for binaries is the 
Deep Ecliptic Survey {Millis et al., 2002; Elliot et al., 
2005). Kern and Elliot (2006) searched 634 unique objects 
from the DES survey and identified 1 . These observations 
were made with 4m telescopes at CTIO and KPNO utiliz- 
ing wide field Mosaic cameras {Muller et al, 1998) with 
0.5 arcsec pixels. Median seeing for the entire data set was 
1.65 arcsec {Kern, 2005) and effectively sets the detection 
limit on angular separation. Magnitude limits in the broad 
VR filter for a well-separated secondary vary, depending 
on seeing, from VR = 23 to VR = 24. 

Follow-up observations of 212 DES objects made with 
the Magellan telescopes to improve astrometry were also 
searched for undetected binaries {Kern and Elliot, 2006). 
Most observations were made with the MagIC camera {Osip 
et al, 2004) which has a pixel scale of 0.069 arcsec. Kern 
(2005) reports the median seeing for these observations was 
0.7 arcsec with a magnitude limit similar to the DES survey. 
Of the 212 objects observed with Magellan, 3 were found to 
be binaries {Osip et al. , 2003, Kern and Elliot, 2005, 2006). 

The Keck telescope survey reported by Schaller and 
Brown (2003) observed 150 objects and found no new bina- 
ries. The observational Umits for this survey have not been 
pubUshed; they are probably similar to the DES as the pri- 
mary limiting factor is seeing. 

The target hsts for the three large ground-based surveys 
have not been pubUshed and it is unclear how much over- 
lap there may be. However, even duplicate observations 
of an individual target can be useful since some TNBs are 
known to have significantly eccentric or edge-on orbits and 
are, therefore, variable in their detectability. The one firm 
conclusion, however, that can be reached from these data 
is that binaries separated sufficiently for detection with un- 
corrected ground-based observations are uncommon, occur- 
ring around 1-2% of TNOs at most. 



The most productive tool for finding TNBs is the HST 
which has found 41 of the 52 companions listed in Table 1. 
The combination of high angular resolution, high sensitiv- 
ity, and stable point- spread-function (PSF) make it ideally 
matched to the requirements for finding and studying TNBs. 
The first conscious search for TNBs using HST was carried 
out in August 2000-January 2001 in a program that looked 
at just 2 TNOs and found no companions. Two other small 
programs executed between October 1997-September 1998 
observed 8 TNOs with the potential to identify a binary, had 
there been one among the objects observed. 

The first moderate-scale program to search for binaries 
used STIS in imaging mode to search for binaries around 
25 TNOs from August 2001 -August 2002. This program 
found 2 binaries, both relatively faint companions to the res- 
onant TNOs (47171) 1999 TCse and (26308) 1998 SM165 
(Trujillo and Brown, 2002; Brown and Trujillo, 2002). The 
STIS imaging mode has a pixel scale of 50 milliarcsec mak- 
ing it possible to directly resolve objects separated by ap- 
proximately 100 milliarcsec or more. In principle, PSF 
analysis can detect binaries at significantly smaller sepa- 
rations in HST data because of the stabiUty of the PSF 
{Stephens and Noll, 2006). The STIS data were obtained 
without moving the telescope to track the target's motion. 
The TNOs observed this way drift measurably during expo- 
sures, complicating the PSF analysis for these data. 

From July 2001-June 2002 75 separate TNOs were ob- 
served with WFPC2 in a program designed to obtain V, R, 
and 1 band colors {Stephens et ai, 2003); 3 of these were 
found to be binary. To achieve better sensitivity and because 
of the relatively large uncertainties in TNO orbits at the 
time, the targets were observed with the WF camera with 
100 milharcsec pixels. The sensitivity to faint companions 
for these data has not been fully quantified and exposure 
times were varied depending on the anticipated brightness 
of the TNO so that, in any case, the hmits vary. Typical pho- 
tometric uncertainties ranged from 3-8% for V magnitudes 
that ranged from 21.9 to 25.1 with a median magnitude of 
23.6. 

NICMOS was used to observe 82 separate TNOs from 
August 2002 through June 2003. Observations were made 
with the NIC2 camera at a 75 milliarcsec pixel scale. Two 
broad filters, the FllOW and F160W, approximating the J 
and H band filters, were used in the observations. A to- 
tal of 9 new binary systems were identified from this data 
set, 3 of which were resolved and visible in the unpro- 
cessed data. The other 6 binaries were identified from a 
process of PSF fitting that enabled a significant increase 
in detectivity beyond the usual Nyquist Umit with a sep- 
aration/relative brightness/total brightness function that is 
complex and modeled with simulated binaries. Several of 
the binaries identified in this way have been subsequently 
resolved with the higher resolution HRC, verifying the ac- 
curacy of the analysis of the NICMOS data {Stephens and 
Noll, 2006). 

From July 2005 through January 2007, HST's HRC was 
used to look at more than 100 TNOs. This program identi- 
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fied a significant number of new binaries {Noll et al, 2006 
a,b,c,d,e,f). The pixels of the HRC, while operational, were 
a distorted 25 x 28 milliarcsec quadrilateral. Observations 
with the Clear filter were able to reach a faint limiting mag- 
nitude of at least V=27, significantly deeper than most other 
HST observations of TNOs (Noll et al., 2006h). 

The wide variety of instruments and inherent observing 
limitations for each makes analysis from a concatenation 
of these data extremely problematic. Indeed, even obser- 
vations taken with the same instrument vary in their lim- 
its based on exposure time, focus, seeing, and other sys- 
tematics. This variability is reduced for space-based obser- 
vations, but not entirely eliminated. Additionally, the de- 
tection of close pairs depends on the ability to model and 
subtract the PSF of the primary with a resulting spatially- 
dependent detection limit. 

2.3 Binary Frequency 

Does the fraction of TNBs vary with any of the observ- 
able properties of TNOs? In the search for correlations, the 
most promising quantities are those that can be associated 
with formation or survival. 

Dynamical class is one such quantity because objects in 
different dynamical classes may have had different origins 
and dynamical histories. On one extreme, objects on un- 
stable, non-resonant, planet-crossing orbits (the Centaurs 
and similar objects) have many close encounters with gi- 
ant planets during their lifetimes. These close encounters 
can potentially disrupt weakly bound binaries (Noll et al., 
2006g; Petit and Mousis, 2004; see also §3.3). At the 
other extreme, objects in the classical disk may have per- 
sisted largely undisturbed from the protoplanetary disk and 
thus be a more congenial environment for survival of bina- 
ries. Stephens and Noll (2006) have shown that, indeed, 
the binary frequency in the cold classical disk is signifi- 
cantly higher than for other dynamical classes at the obser- 
vational limits attainable by the NIC2 camera. They found 
that 22± of classical TNOs with inclinations less than 

5 

5° had detectable binaries at magnitude-dependent separa- 
tions typically > 60 milliarcsec. The rate of binaries for 
all other dynamical classes combined was 5.5±2% for the 
same limits. The number of objects and binaries in this 
dataset precluded any meaningful conclusions on the binary 
frequency in the various dynamically "hot" populations. In- 
cluding more recent HRC observations strengthens these 
conclusions, as shown in Figure 2. 

The size of the primary is another parameter that could, 
potentially, be correlated with binary status. To the extent 
that size correlates with dynamical class, the two sorting 
criteria can be confused for one another. There are, in fact, 
differences in size distributions as a function of dynami- 
cal class with the resonant objects, scattered disk, and high 
inclination classical objects all having significantly larger 
upper size limits than the cold classical population (Levison 
and Stern, 2001; Bernstein et al, 2004;). The published 
binary searches to date cover an insufficient range in size 



(as measured by Hv) to be able to reach any strong conclu- 
sions with regard to binary frequency as a function of size. 
Brown et al. (2006a) have noted an apparently higher frac- 
tion of bound systems among the largest TNOs based on 
the satellites known for Pluto, Eris, and 2003 ELgi. With 
the size bins chosen in this work (the four largest TNOs 
versus smaller TNOs) the difference in binary frequency 
is mathematically significant. The subsequent detection of 
small satellites around four more large TNOs {Brown and 
Suer, 2007) strengthens this apparent trend. In considering 
whether small satellites exist around smaller TNOs, it is im- 
portant to consider that companions as faint as those of the 
largest TNOs would only be detectable for widely separated 
companions in a subset of the most recent deep observa- 
tions (Figure 3 Noll et al, 2006g) and would not have been 
detectable in most earlier surveys. At the same fractional 
Hill radius, small satellites of 100-km class TNOs would 
be extremely difficult to detect in any existing observations. 
Caveats aside, however, it seems reasonable to hypothesize 
that the small satellites of the largest TNOs may be colli- 
sional in origin {Brown et al, 2006a; Stern et al, 2006) 
while the nearly equal sized binaries of smaller TNOs may 
form from dynamical capture and thus, there may be real 
differences in the frequencies of these two types of bound 
systems. We discuss these two modes of formation in more 
detail in §4. 

An obvious, but sometimes neglected, quahfier that must 
accompany any description of binary frequency is the limit 
in magnitude and separation imposed by the observational 
method used to obtain the data. This is especially true be- 
cause, as discussed in more detail below and as shown in 
Figure 4, the number of binaries detectable in a given sam- 
ple appears to be a strong function of separation {Kern and 
Elliot, 2006). The variation of binary frequency with dy- 
namical class and/or size means that statements about bi- 
nary frequency must be further qualified by a description 
of the sample. Thus, it is impossible to state a unique 
"frequency of transneptunian binaries". The need for a 
more nuanced description of binary frequency among the 
transneptunian populations is both a challenge and an op- 
portunity waiting to be exploited. 

3. PHYSICAL PARAMETERS 

3.1 Relative Sizes of Binary Components 

The relative sizes of the primary and secondary com- 
ponents is an important physical parameter that is usually 
available from a single set of observations (with the impor- 
tant caveat of possible non-spherical shapes and lightcurve 
variations as we discuss in §3.6). The observed magnitude 
difference between the two components of a binary, A^ag, 
can be used to constrain the ratios of their radii R\ / R2 and 
surface areas A1/A2 according to 

i?2 \ A2 V Pi 
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Table 1 
Transneptunian Binaries 
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r9] 


2000 CQ114 


c 


46, 


.230 


0, 


.110 


2.7 


0.178(5) 


0.4(2) 


6.6 


no] 


2000 CF105 


c 


43, 


.881 


0, 


.037 


0.5 


0.78(3) 


0.7(2) 


6.9 


rii] 


1999 OJ4 


c 


38, 


.067 


0, 


.023 


4.0 


0.097(4) 


0.16(9) 


7.0 


r2] 


2001 FL186 


c 


44, 


.178 


0, 


.077 


3.6 


0.065(14) 


0.8(6) 


7.0 


rti 


2003 UN284 


c 


42, 


.453 


0, 


.010 


3.1 


2.0(1) 


0.6(2) 


7.4 


ri2] 


2001 QW322 


c 


44, 


.067 





.027 


4.8 


4 


0.0(1) 


7.8 


ri3] 


1999 RT214 


c 


42, 


.711 


0, 


.052 


2.6 


0.107(4) 


0.81(9) 


7.8 


[14] 


2003 TJss 


c 


44. 


.575 


0, 


.089 


1.0 


0.119(2) 


0.50(7) 


7.8 


[tl 


Scattered 






















(136199) Eris 


s 


67. 


.695 


0, 


.441 


44.2 


0.53(1) 


4.43(5) 


-1.2 


[15] 


(136108) 2003 ELei 


s 


43. 


.316 


0, 


.190 


28.2 


0.63(2) 


3.1(1) 


0.2 


[16] 
















0.52(3) 


4.6(4) 




[17] 


(55637) 2002 UX25 


s 


42, 


.551 





.141 


19.5 


0.164(3) 


2.5(2) 


3.6 


[1] 


(120347) 2004 SBgo 


s 


42, 


.032 


0, 


.104 


24.0 


0.107(3) 


2.2(1) 


4.4 


[18] 


(48639) 1995 TLg 


X 


52, 


.267 


0, 


.235 


0.2 


0.01(1) 


1.7 


5.4 


r2] 


2001 QC29S 


s 


46, 


.222 





.123 


30.6 


0.130(7) 


0.58(3) 


6.1 


[19] 


2004 PB108 


s 


44, 


.791 


0, 


.096 


20.3 


0.172( 3) 


1.2(1) 


6.3 


[20] 


(65489) Ceto/Phorcys 





102.876 


0, 


.821 


22.3 


0.085(2) 


0.6(1) 


6.3 


[21] 


2002 GZ31 


X 


50, 


.227 


0, 


.237 


1.1 


0.070(9) 


1.0(2) 


6.5 


[22] 


(60458) 2000 CMim 


s 


59. 


.838 


0, 


.407 


19.7 


0.074(6) 


0.57(7) 


7.0 


[23] 


(42355) Typhon/Echidna 





38. 


.112 


0, 


.540 


2.4 


0.109(2) 


1.47(4) 


7.2 


[24] 


Resonant 






















(134340) Pluto/Charon 


3:2 


39. 


.482 


0, 


.248 


17.1 


0.9 


2-3 


-1.0 


[25] 


Nix 














1.85(4) 


9.26(2) 




[26] 


Hydra 














2.09(4) 


8.65(2) 




[26] 


(90482) Orcus 


3:2 


39. 


.386 


0, 


.220 


20.6 


0.256(2) 


2.5 


2.3 


[1] 


2003 AZ84 


3:2 


39. 


.414 


0, 


.181 


13.6 


0.22(1) 


5.0(3) 


3.9 


[1] 


(47171) 1999 TC36 


3:2 


39. 


,270 


0, 


.222 


8.4 


0.367(4) 


2.21(1) 


4.9 


[27] 


(82075) 2000 YW134 


8:3 


57. 


.779 


0, 


.287 


19.8 


0.06(1) 


1.3 


5.0 


[2] 


(119979) 2002 WC19 


2:1 


47. 


,625 


0, 


.260 


9.2 


0.090(8) 


2.5(4) 


5.1 


[28] 
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where pi and p2 are the albedos of the TNB components. It 
is usual, but not necessary, to assume that both components 
have the same albedo. As we note in §3.6, the similarity of 
the colors of TNB components suggests connmon surface 
materials with similar albedos may be the norm for TNBs. 
However, we also note that in the one instance where sep- 
arate albedos have been measured, the Pluto/Charon sys- 
tem, they are not identical. Given the large range of albedos 
of TNOs of all sizes {e.g. Grundy et ai, 2005; chapter by 
Stansberry et ai), it is reasonable to keep this customary 
simplification in mind. 

The relative sizes of TNBs found so far is heavily skewed 
to nearly-equal sized systems as can be seen in Table 1 and 
Figure 3. The prevalence of nearly equal-sized systems is 
a unique feature of TNBs compared to binaries in the Main 
Belt or Near Earth populations (e.g. Richardson and Walsh, 
2006; Noll, 2006). To some extent, this conclusion is lim- 
ited by observational bias, but, deep surveys with the HST's 
HRC show an apparent lack of asymmetric binaries {Noll et 

ah, 2006h; Figure 3). A preference for similar sized compo- 
nents is a natural outcome of dynamical capture models for 
the formation of binaries {Astakhov et ai, 2005) and may 
explain the observed distribution of relative sizes after ac- 
counting for the small number of large object satellites that 
appear to have been formed from collisions {e.g. Stern et 

ai, 2006). 

3.2 Orbit Determination 

Determining the mutual orbit of a spatially-resolved pair 
of gravitationally bound objects is a classic problem of 
celestial mechanics. Solutions are reviewed in numerous 
textbooks {e.g. Aitken, 1964; Heintz, 1978; Smart, 1980; 
Danby, 1992). 

The Keplerian orbit of a pair of bound point masses can 
be described by seven independent quantities. There is con- 
siderable flexibility in the choice of the seven. One possible 
set is the system mass plus the three dimensional relative 
position and velocity vectors at a specific time. This form is 
generally preferred for specifying inputs to orbital integra- 
tion routines. A second set of elements is more frequently 
used to specify binary orbits in the scientific hterature: pe- 
riod T, semimajor axis a, eccentricity e, inclination i, mean 
longitude at a specific time e, longitude of the ascending 
node n, and longitude of periapsis w. Other sets of ele- 
ments are possible and are occasionally encountered in the 
literature. 

Each observation of a resolved binary pair provides two 
constraints, the separation and position angle, or equiva- 
lently, the relative positions in right ascension and dec- 
hnation. To constrain seven orbital elements requires at 
least four observations, in order to have more constraints 
than unknowns. In practice, four observations are often not 
enough to uniquely determine the seven unknowns, because 
the constraints from the separate observations are not nec- 
essarily independent of one another. For example, four ob- 
servations that happen to sample the same orbital longitude 



do little to constrain the shape and orientation of the orbit. 
Many more than four observations may be needed, if they 
are poorly timed. 

Instruments capable of doing relative astrometry on the 
extremely close-spaced and faint components of TNBs are 
a scarce and valuable resource. To make the most efficient 
use of these facilities requires strategic timing of follow- 
up observations to maximize the additional constraints pro- 
vided by each successive observation. Several groups have 
applied Monte Carlo techniques to this task {e.g., Margot et 
al, 2004, 2005; Hestrojfer et at. 2005; Grundy et al. 2007). 
The general approach is to generate random orbits consis- 
tent with the existing set of observations which do not yet 
uniquely determine the binary orbit. The collection of or- 
bits produced by this exercise is used to map out regions 
of orbital element space consistent with the constraints al- 
ready available, and to identify times when follow-up ob- 
servations would be most effective for collapsing the cloud 
of possible orbits. 

An additional complication involves near-equal bright- 
ness binaries, which are not uncommon among the TNBs. 
When the "primary" and "secondary" have similar bright- 
ness, or large amplitude lightcurve variations, they become 
difficult to distinguish from one another, leading to uncer- 
tainties which require additional observations to resolve. 
The binary system 2003 QYgo provides a recent example 
of this situation {Kern and Elliot, 2006). 

Various non-linear least squares techniques such as the 
Levenburg-Marquardt and AMOEBA algorithms {Press et 
al, 1992) can be used to fit a set of orbital elements to the 
observational data by iteratively minimizing the residuals 
(as measured by the statistic) between the astrometric 
data and the positions computed from the orbital elements. 
In situations where observational data are particularly abun- 
dant, linear least squares fitting techniques may be appli- 
cable. Knowledge of the uncertainties in the fitted orbital 
elements is just as important as knowledge of the elements 
themselves. The Monte Carlo techniques mentioned earlier 
can be used to directly investigate uncertainties in orbital 
elements {e.g., Virtanen etal., 2001, 2003; Hestrojfer et al., 
2005). Another approach involves varying each parameter 
around its best fit value and allowing the remaining param- 
eters to readjust themselves to compensate, resulting in a 
new fit with different y^. This approach can be used to map 
out a seven-dimensional space, within which contours of 
1 fj or 3 cr likeUhood can be traced {Lampton, 1976; Press 
etal, 1992). 

For TNBs, relative motion between the Earth and the or- 
biting pair leads to changing viewing geometry over time. 
The motion has two components: the comparatively rapid 
motion of the Earth around the sun and the much slower 
motion of the TNB along its heliocentric orbit. The re- 
sulting changes in viewing geometry complicate the orbit- 
fitting process compared with procedures developed for bi- 
nary stars which assume fixed geometry between the ob- 
server and the binary system. However, there are advan- 
tages to observing a binary from different angles. Relative 
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Table 1 — Continued 



dyn. heliocentric orbit elements separation 
object class a© (AU) e© i© (°) so (arcsec) Amag ref 



(26308) 1998 SMies 


2:1 


47.501 


0.370 


13.5 


0.205(1) 


2.6(3) 


5.8 


[29] 


2003 QWiii 


7:4 


43.659 


0.111 


2.7 


0.321(3) 


1.47(8) 


6.2 


[30] 


2000 QL251 


2:1 


47.650 


0.216 


3.7 


0.25(6) 


0.05(5) 


6.3 


[31] 


(60621) 2000 FEs 


5:2 


55.633 


0.404 


5.9 


0.044(3) 


0.6(3) 


6.7 


[20] 


(139775) 2002 QG298 


3:2 


39.298 


0.192 


6.5 


contact? 


N/A 


7.0 


[32] 



Note. — Objects are sorted into three dynamical groups. Classical (both hot, H, and cold C), Scattered 
(includes Scattered-near, S, Scattered-extended, X, and Centaurs, 0) and Resonant. Objects in each grouping 
are sorted by absolute magnitude, Hv- Uncertainties in the last digit(s) of measured quantities appear in 
parentheses. Table entries in italics indicate quantities that have been published without error estimates or 
that have been computed by the authors from estimated quantities. The Hy colimm Usts the combined 
absolute magnitude of the system as tabulated by the Minor Planet Center (MFC). * The brightness of Logos 
varies significantly as described in §3.6. References: [1] Brown and Suer (2007); [2] Stephens and Noll 
(2006); [3] Elliot et al. (2001); [4] Noll et al. (2003); [5] Veillet et al. (2001); [6] Kern and Elliot (2005); 
[7] Elliot, Kern, and Clancy (2003); [8] Noll et al. (2007a); [9] Noll et al. (2002a); [10] Stephens, Noll, and 
Grundy (2004); [11] Noll et al. (2002b); [12] Millis and Clancy (2003); [13] Kavelaars et al. (2001); [14] 
Noll et al. (2006f); [15] Brown (2005a); [16] Brown et al. (2005b); [17] Brown (2005b); [18] Noll et al. 
(2006e); [19] Noll et al (2002c); [20] Noll et al (2007d); [21] Noll et al. (2006g), Grundy et al (2006); 
[22] Noll et al (2007c); [23] Noll et al. (2006b); [24] Noll et al. (2006a); [25] Smith et al. (1978), Christy 
and Harrington (1978); [26] Weaver et al (2005); [27] Trujillo and Brown (2002); [28] Noll et al. (2007b); 
[29] Brown and Trujillo (2002); [30] Noll et al. (2006c); [31] Noll et al. (2006d); [32] Sheppard and Jewitt 
(2004);[t] unpubUshed as of 28 February 2007. 
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Fig. 2. — Classical TNOs observed by HST are plotted with their absolute magnitude, Hv, on the horizontal axis and their inclination 
to the ecliptic, i, plotted vertically. Objects observed with the HRC ( 30 milliarcsec resolution) are shown as large circles; less sensitive 
observations made with NIC2 ( 75 milliarcsec resolution) are shown as smaller circles. Binaries are shown as filled circles. The 
dotted Une at 4.6° is the boundary between Hot and Cold Classical populations proposed by Elliot et al. (2005). The extremely strong 
preference for low inclination binaries in this sample is evident. 
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Fig. 3. — In this figure the y-axis is the observed magnitude of the primary component of observed TNBs. The x-axis shows the range 
of observed and detectable magnitude differences, Amag, for secondaries under several different observing circumstances. Binaries 
detected with HST's HRC are shown as triangles, binaries found with other instruments are shown as diamonds. The stippled area shows 
the 7-sigma detection limit for companions that are sufficiently separated that their local background is dominated by sky noise and 
dark current. Objects that are observed by the HRC and are background dominated are shown as large triangles. The background and 
detection limits for the objects shown by the small triangles are limited to a varying degree by the PSF of the primary. The dashed line 
defines an approximate empirical detection limit for objects separated by 3-pixels from their primary, i.e. 75 milliarcsec. The backgroimd 
at this separation is dominated by the PSF of the primary to a degree that varies as a function of the primary. Both detection limits are 
dependent on the details of the observations and therefore do not apply to the observations from other instruments (diamonds) which are 
less sensitive. The clustering of binaries at Amag < 1 appears to be an intrinsic property of TNBs and not an observational bias. (Figure 
from Noll et ai, 2006h) 



astrometry is otherwise not able to distinguish between the 
actual orbit and its mirror reflection through the sky plane. 
For TNBs, the ambiguity between these two "mirror" orbits 
can be broken by observations spanning as Uttle as a few 
years {e.g. Hestrojfer and Vachier, 2006), but this has only 
actually been done for a small minority of orbits to date. 

For larger and closer binaries, departures from spherical 
symmetry of the primary's mass distribution (parametrized 
by the gravitational harmonic coefficients, primarily J2) can 
exert torques on the secondary leading to secular preces- 
sion of the line of nodes and/or of the line of apsides {e.g., 
Brouwerand Clemence, 1961). Observations of this secular 
evolution can be used to measure J2 as has been done for a 
binary asteroid {Marchis et ai, 2005), but this has not yet 
been done for any TNBs except Pluto (and that measure- 
ment by Descamps, 2005 neglected possible perturbations 
on the orbit of Charon by Pluto's smaller satellites). Com- 
parable effects can be anticipated for the orbits of the more 
distant satellites of close pairs, such as the orbits of Hydra 
and Nix around the Pluto/Charon binary {Buie et al. 2006; 
Lee and Peale, 2006). 



3.3 Semimajor Axis Distribution 

The separation of components in binaries is a fundamen- 
tal datum that can be diagnostic in comparing systems and 



in understanding general processes affecting populations. 
Semimajor axes are known for only the subset of TNBs with 
measured orbits (Table 2). Unfortunately, from a statistical 
perspective, the available data are too few and are heavily 
biased against small semimajor axes by limitations of both 
discovery and follow-up observations. 

As a proxy for the semimajor axis, we can use the sepa- 
ration at discovery, ,so, listed in Table 1. In any individual 
case this is only an approximate substitute for a measure- 
ment of the semimajor axis from an orbit because of the 
unknown orientation of the orbit plane with respect to the 
line of sight and the unknown eccentricity of the binary or- 
bit. Observational biases tend to favor detection of objects 
close to their widest separation from the primary indepen- 
dent of the orientation of the orbit plane. For the 13 systems 
in Table 2 the ratio of the semimajor axis to the distribu- 
tion of separation at discovery, a/so, ranges from 0.5 to 
2. The median value of u/sq is 0.96. Thus, for a moder- 
ate size ensemble, the distribution of separations at discov- 
ery provides an excellent statistical approximation for the 
semimajor axis distribution. The separations at discovery 
from Table 1 are plotted in Figure 4. This log plot shows 
a stronger than exponential increase at smaller separations. 
The trend of increasing binary frequency at small separation 
is robust because observational bias will decrease the num- 
ber of detectable companions near the limits of instrument 
resolution. Interestingly, this trend is in qualitative agree- 
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Table 2 

Physical Properties of Transneptunian Binaries 



object 


a 

(km) 


e 


T 

(days) 


M 
(10^* kg) 


p\ 


p 

(g cm"'^) 


J/J' 


refs 


(136199) Eris 


36,000 





14 


16,400 


0.86(7) (V) 


2.26(25) 


0.16 


[1] 


(134340) Pluto/Charon 








14,570(9) 






0.40 


[2] 


Pluto 








13,050(620) 


0.51-0.71 (V) 


2.03(6) 




[2] 


Charon 


19,571(4) 


0.00000(7) 


6.387230(1) 


1,521(65) 


0.38 (V) 


1.65(6) 




[4] 


Nix 


48,670(120) 


0.002(2) 


24.8562(13) 


0.1-2.7 


0.01-0.35 


2.0 




[4] 


Hydra 


64,780(90) 


0.005(1) 


38.2065(14) 


0.2^.9 


0.01-0.35 


2.0 




[4] 


(136108) 2003 ELei 


49,500(400) 


0.050(3) 


49.12(3) 


4,210(100) 


0.7(1) iv) 


2.9(4) 


0.53 


[3] 


(47171) 1999 TC36 


7,720(460) 


0.22(2) 


50.4(5) 


14.4(2.5) 


0.08(3) (v) 


0.5(3/2) 


0.31 


[5] 


2001 QC298 


3,690(70) 


0.34(1) 


19.2(2) 


10.8(7) 


0.08 (Vfm) 


1.0 


1.16 


[6] 


(26308) 1998 SMies 


11,310(110) 


0.47(1) 


130(1) 


6.78(24) 


0.08(3/2) (V) 


0.7(3/2) 


0.56 


[7] 


(65489) Ceto/Phorcys 


1840(50) 


<0.014 


9.557(8) 


5.4(4) 


0.08(2) (V) 


1.4(6/3) 


0.89 


[8] 


(66652) 1999 RZ253 


4660(170) 


0.46(1) 


46.263(6/74) 


3.7(4) 


0.17 (R) 


1.0 


1.56 


[9] 


1998 WW31 


22,300(800) 


0.82(5) 


574(10) 


2.7(4) 


0.054 (R) 


1.0 


2.22 


[10] 


(88611) 2001 QT297 


27,880(150) 


0.241(2) 


825(1) 


2.51(5) 


0.13 (r') 


1.0 


3.41 


[11] 


(42355) Typhon/Echidna 


1830(30) 


0.53(2) 


18.971(1) 


0.96(5) 


0.05(V) 


0.47(18/10) 


2.13 


[12] 


(58534) Logos/Zoe 


8,010(80) 


0.45(3) 


312(3) 


0.42(2) 


0.37(4) (R) 


1.0 


2.65 


[13] 


2003 QY90 


7,000-13,000 


0.44-0.93 


306-321 


0.3-1.7 


0.19-0.35 (t'} 


1.0 


2.14 


[14] 



Note. — Objects are listed in order of decreasing system mass. Uncertainties in the last digit(s) of measured quantities appear in parentheses. 
Table entries in italics indicate quantities that have been published without error estimates or that have been computed by the authors from 
estimated quantities. The components of the Pluto system are listed separately because they have had independent determinations of radius, 
mass, albedo, and density that are not yet available for other binary or multiple systems. The wavebands used to determine geometric albedos 
are shown in parentheses {v is for references citing "visual" albedo, Veoe is the HRC's F606W filter). Densities of 1.0 g cm^'^ are assumed 
for objects without an independent size measurement. The corresponding geometric albedos are for this assumed unit density. References: 
[1] Brown et al. (2006a,b); [2] Buie et al. (1997), Buie et al. (2006), Rabinowitz et al. (2006); [3] Brown et al. (2005); [4] Buie and Grundy 
(2000), Buie et al. (2006), Lee and Peale (2006); [5] Stansberry et al. (2006); [6] Margot et al. (2004); [7] Margot et al. (2004), Spencer et al. 
(2006); [8] Grundy et al. (2007); [9] Noll et al (2004a); [10] Veillet et al. (2002); [11] Osip et al. (2003), Kern (2005); [12] Grundy et al, in 
press; [13] Noll et al. (2004b), Margot et al. (2004); [14] Kern et al., in preparation 
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Fig. 4. — The number of objects with a separation greater than or equal to a given separation is shown. The number of objects at small 
separation increases faster than an exponential as can be seen in this logarithmic plot. Observational bias favors detection of widely 
separated binaries suggesting that the underlying distribution may rise even more rapidly with decreased separation. 



ment with the binary formation model discussed in more 
detail in §4.5. 

It is natural to wonder if the trend of binary frequency 
continues at even smaller separations or if there is some 
eventual cutoff. The separation of TNB components in 
terms of "typical" values for the distance to Earth in AU, 
A®, and an angular separation in arcsec, 6, can be ex- 
pressed as 

or nearly 60 primary radii for a 100 km diameter TNO. This 
leaves a significant amount of separation phase space where 
a sizable population of stable binaries could exist undiscov- 
ered. 

At the extreme of small separation are bilobed objects 
and contact binaries. Sheppard and Jewitt, (2004) proposed 
that the TNO 2001 QG298 may be a very close, possibly 
contact binary based on its large amplitude lightcurve (1.14 
mag) and its relatively slow rotation (13.77 hr). The rota- 
tion is too slow and the ampUtude too large to explain the 
lightcurve as a fluid deformation into a Jacobi ellipsoid as 
is postulated for some other TNOs (Sheppard and Jewitt, 
2002; Jewitt and Sheppard, 2002). That leaves only albedo 
variation, strength-dominated shape, or a close binary to ex- 
plain the lightcurve of 2001 QG298. Sheppard and Jewitt, 
(2004) conclude that a close binary is the most likely expla- 
nation for the observed lightcurve. They further extrapolate 
from this single object to suggest that as many as 10-15% 
of TNOs could be close binaries. While unproven, a 10- 



15% fraction of TNOs as contact binaries is not incompat- 
ible with the frequency of wider binaries and the trend of 
frequency as a function of separation. 

The natural dimension for scaling the separations of 
TNBs is the Hill radius, i?H, given by the equation 

i?H ^ a0(Mi/3MQ)V3 , (3) 

where Mi is the mass of the largest component of the bi- 
nary, Mq is the mass of the Sun, and a© is the semimajor 
axis of the primary's heliocentric orbit. SatelUte orbits well 
inside the primary's Hill sphere are generally stable with re- 
spect to external perturbations. However, the instantaneous 
Hill radius, calculated by replacing a© above with the in- 
stantaneous heliocentric distance, can deviate significantly 
from the average Hill radius for TNBs with large heliocen- 
tric eccentricities such as (65489) Ceto/Phorcys. This result 
must be taken into account when considering, for example, 
the impact of gravitational perturbations from major plan- 
ets. The Hill radius scales with the radius of the primary 
and, thus, the large size range of the primary objects in Ta- 
ble 2 implies a similarly large range in Hill radii measured 
in absolute units. In dimensionless units, a typical Hill ra- 
dius for a TNB is on the order of 7000 times the radius of the 
primary. Objects with measured orbits have semimajor axes 
that range from -^0.1% to ~8% of the Hill radius, {Noll et 
al, 2004b; Kern, 2005) well within the stable portion of the 
Hill sphere. Whether this represents an intrinsic outcome 
of formation {Astakhov et al, 2005), or is a signature of 
a surviving remnant population where more weakly bound 
systems have been disrupted while more tightly bound sys- 
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tems have become even tighter in the wake of encounters 
with third bodies {Petit and Mousis, 2004), or some combi- 
nation of the two, remains to be determined. 

3.4 Mass, Albebo and Density 

A particularly valuable piece of information that can be 
derived from the mutual orbit of a binary system is the total 
mass, Mgys, of the system, according to the equation 

Msys = , (4) 

where a is the semimajor axis, G is the gravitational con- 
stant, and T is the orbital period. Knowledge of a tends to 
be limited by the spatial resolution of the telescope, while 
knowledge of T is limited by the timespan over which 
observations are carried out (modulo the binary orbit pe- 
riod). It is possible to extend the timespan of observa- 
tions, whereas the spatial resolution cannot generally be im- 
proved. Thus, typically, T is determined to much higher 
fractional precision than a, and the uncertainty in A/gys is 
dominated by the uncertainty in a. Mgys can often be cal- 
culated before all 7 elements of the binary orbit are fully 
determined, since the elements T, a, and e are relatively in- 
sensitive to the ambiguity between orbits mirrored through 
the instantaneous sky plane (see §3.1). 

For a system with a known mass, it is possible to derive 
other parameters which offer potentially valuable composi- 
tional constraints. For instance, if one assumes a bulk den- 
sity p, the bulk volume of the system Vgys can be computed 
as Vsys = Msyjp. How the volume and mass is partitioned 
between the two components remains unknown. Assuming 
the components share the same albedo, the individual radii 
of the primary and secondary can be obtained from 

„ ( SVsys 

^ \4Tr{l - lO-o.6A„„,) J 

and equation ^ simplifies to give i?2 = i?ilO^°'^'^"'*«. 

An effective radius Rcb, equal to the radius of a sphere 
having the same total surface area as the binary system can 
be computed as 

i?cff = ^Rl + Rl. (6) 

Combining i?off and the absolute magnitude of the sys- 
tem Hx, one obtains the geometric albedo p\ 

where C\ is a wavelength-dependent constant (Bowell, 
1989; Harris, 1998). For observations in the V band 
Cy = 664.5 km. This approach has been used to estimate 
albedos for a number of TNBs by assuming their bulk den- 
sities must lie within a plausible range, typically taken to be 
0.5 to 2 g cm"3 (e.g., Noll et al, 2004a, 2004b; Grundy et 



ai, 2005; Margot et al, 2005). These efforts demonstrate 
that the TNBs have very diverse albedos, but those albe- 
dos are not clearly correlated with size, color, or dynamical 
class. The calculation could also be turned around such that 
an assumed range of albedos leads to a range of densities. 

When the sizes of the components of a binary system 
can be obtained from an independent observation, that in- 
formation can be combined with the system mass to obtain 
the bulk density, providing a fundamental constraint on bulk 
composition and interior structure. Sizes of TNOs are ex- 
tremely difficult to obtain, although a variety of methods 
can be used, ranging from direct observation {e.g.. Brown 
et al., 2004) to mutual events and stellar occultations {e.g., 
Gulbis et al., 2006). For rotationally deformed bodies it 
is possible to constrain the density directly from the ob- 
served lightcurve assuming the object is able to respond as 
a "fluid" rubble pile {Jewitt and Sheppard, 2004; Takahashi 
and Ip, 2004; chapter by Sheppard et al.). Spitzer Space 
Telescope observations of thermal emission have recently 
led to a number of TNO size estimates {e.g., Cruikshank 
et al, 2006; chapter by Stansberry et al.). Unfortunately, 
many of the known binaries are too small and distant to be 
detected at thermal infrared wavelengths by Spitzer or di- 
rectly resolved by HST. 

Systems with density estimates include three large 
TNBs: Pluto and Charon, with p = 2.0±0.06 and p = 
1.65±0.06 g cm-3, respectively {Buie et al, 2006), 2003 
ELgi with p = 3.0±0.4 g cm^'' {Rabinowitz et al, 2006), 
and Eris with p = 2.26±0.25 g cm^^ (^Brown, 2006). The 
relatively high densities of the large TNOs are indicative 
of substantial amounts of rocky and/or carbonaceous ma- 
terial in the interiors of these objects, quite unhke their 
ice-dominated surface compositions. 

Four smaller TNBs have recently had their densities de- 
termined from Spitzer radiometric sizes. These are (26308) 
1998 SMi65 with p = 0.70 ± [J ^i g cm"^ {Spencer et al 
2006), (47171) 1999 TC36 with p = 0.5 ± [J;^ g cm-^ 
{Stansberry et al, 2006), (65489) Ceto/Phorcys with p = 
1.38 ± °-32 g cm-3 {Grundy et al, 2007), and (42355) Ty- 
phon/Echidna with p = 0.47 ± [J J* g cm^'^ {Grundy et al 
in preparation). Takahashi and Ip (2004) estimate a density 
of p <0.7 g cm~'^ for 2001 QG298. The very low densities 
of four of these five require little or no rock in their interi- 
ors, and even for pure H2O ice compositions, call for con- 
siderable void space. The higher density of Ceto/Phorcys is 
consistent with a mixture of ice and rock. It is clear from 
these results that considerable diversity exists among den- 
sities of TNOs, but it is not yet known whether densities 
correlate with externally observable characteristics such as 
color, lightcurve amplitude, or dynamical class. 

3.5 Eccentricity Distribution and Tidal Evolution 

The eccentricities of binary orbits known to date span 
the range from values near zero to a high of 0.8 (Table 2), 
with perhaps a clustering in the range 0.3-0.5. With the 
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usual caveats about small number statistics (N^IO), no ob- 
vious correlation between eccentricity and semimajor axis 
is present in the data obtained to date. 

Keplerian two-body motion assumes point masses orbit- 
ing one another. The finite size of real binary components 
allows differential gravitational acceleration between nearer 
and more distant parts of each body to stretch them along 
the line connecting them. If their mutual orbit is eccentric, 
this tidal stretching varies between apoapsis and periapsis, 
leading to periodic flexing. The response of a body to such 
flexing is characterized by the parameter Q, which is a com- 
plex function of interior structure and composition {Goldre- 
ich and Soter, 1966; Farinella et ai, 1979). In a body with 
low Q, tidal flexing creates more frictional heating, which 
dissipates energy. A body with higher Q can flex with less 
energy dissipation. For a body having a rotation state dif- 
ferent from its orbital rotation Q is related to the angular 
lag, S, of the tidal bulge behind the line of centers: Q^^ 
= tan(25). As with orbital eccentricity, this situation pro- 
duces time-variable flexing, and thus frictional heating and 
dissipation of energy. Typical values of Q for rocky plan- 
ets and icy satellites are in the 10 to 500 range (Goldreich 
and Soter, 1966; Farinella et ai, 1979; Dobrovolskis et ai, 
1997). 

The energy dissipated by tidal flexing comes from or- 
bital and/or rotational energy, leading to changes in orbital 
parameters and rotational states over time. Tides raised on 
the primary tend to excite eccentricity, while tides raised on 
the secondary result in damping. The general trend is to- 
ward circular orbits, with both objects' spin vectors aligned 
and spinning at the same rate as their orbital motion. The 
timescale for circularization of the orbit is given by 

j fay 
^"'^ - 63Mi V G(Afi + M2) [R2J 

where a is the orbital semimajor axis. Mi is the mass of 
the primary, and M2, R2, and Q2 are the mass, radius, and 
Q of the secondary (e.g., Goldreich and Soter, 1966). It is 
important to recognize that this formula (a) assumes the sec- 
ondary to have zero rigidity, and (b) assumes that the eccen- 
tricity evolution due to tides raised on the primary is ignor- 
able. Neither of these assumptions may be justified, espe- 
cially in the case of near-equal sized binaries (for a thorough 
discussion, see Goldreich and Soter, 1966). The timescale 
^ is sensitive to the ratio of the semimajor axis to the 
size of the secondary. Larger and/or closer secondaries are 
likely to have their orbits circularized much faster than more 
widely- separated systems. For example, the close binary 
(65489) Ceto/Phorcys has a = 1,840 km. Mi = 3.1x10^^ 
kg, A/2 = 1.7x101** kg, and R2 = 67 km (Grundy et ai, 
2007). For Q = 100 (a generic value for soUd bodies), its 
orbit should circularize on a relatively short timescale of the 
order of ^10^ years. A more widely separated example, 
(26308) 1998 SM165 has a = 11,300 km. Mi = 6.5 xlO^^ 
kg, Af2 = 2.4x101"^ kg, and i?2 = 48 km (Spencer et ai, 
2006), leading to a much longer circularization timescale 



of the order of ^^10 years, consistent with the observa- 
tion that it still has a moderate orbital eccentricity of 0.47 
(Margot, 2004). 

Tidal effects can also synchronize the spin rate of 
the secondary to its orbital period (as in the case of the 
Earth/Moon system) and, on a longer timescale, synchro- 
nize the primary's spin rate as well (as for Pluto/Charon). 
The timescale for spin locking the primary (slowing its spin 
to match the mutual orbital period) is given by 

_ QiRlu;i (MiV ( a V 

where loi is the primary's initial angular rotation rate and 
i?i is its radius (Goldreich and Soter, 1966). The initial 
angular rotation rate is not known, but an upper limit is the 
break-up rotation rate, which is within 50 percent of 3.3 
hours for the 0.5 to 2 g cm"'^ range of densities discussed. 
For (65489) Ceto/Phorcys and (26308) 1998 SM165 with 
primary radii (Ri) of 86 km (Grundy et ai, 2007) and 147 
km (Spencer et al, 2006), we find the spinlock timescale to 
be ^lO"' and ■^lO^ years respectively, slightly faster than 
the circularization timescale. The timescale for de-spinning 
the secondary, Tdospin,2, is given by swapping subscripts 1 
and 2 in 

The general case for both tidal circularization and tidal 
despinning in binaries is complex (e.g. Murray and Der- 
mott, 1999). Binaries in the transneptunian population in- 
clude many where the secondary is of comparable size to 
the primary. For these systems, it is not safe to make the 
common assumption that the tide raised by the secondary 
on the primary is ignorable. Furthermore, there are many 
systems where the binary orbit has moderate to high eccen- 
tricity invalidating the simplifying assumptions possible for 
nearly circular orbits. A full treatment of the tidal dynam- 
ics for the kinds of systems we find in the transneptunian 
population would make an interesting addition to the litera- 
ture. In the meantime, observation is likely to lead the way 
in understanding the dynamics of these systems. 

3.6 Colors and Lightcurves 

The spectral properties of TNOs and their temporal vari- 
ation are fundamental probes of the surfaces of these ob- 
jects (see chapters 7-11, this volume). The colors of TNOs 
have long been known to be highly variable (Jewitt andLuu, 
1998) and some correlations of color with other physical or 
dynamical properties have been claimed (e.g. Peixinho et 
ai, 2004). A natural question is whether this variability 
can be used to constrain either the origin of TNBs, the ori- 
gin of color diversity, or both. For example, one can ask 
whether TNB components are similarly or differently col- 
ored. Because TNBs are thought to be primordial, differ- 
ences in color could be either due to mixing of different- 
composition populations in the protoplanetary disk before 
the bound systems were formed or different collisional and 
evolutionary histories of components after they are bound. 
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A handful of TNBs have reported single-epoch resolved 
color measurements. Some of these objects are solar- 
colored [2000 CF105, (58534) Logos/Zoe, (47171) 1999 
TC36, (66652) 1999 RZ253, (88611) 2OOIQT297], while 
others are more gray [2OOIQC298, (65489) Ceto/Phorcys]. 
However, so far, the components have colors that are consis- 
tent with each other within the uncertainties of the measure- 
ments, 0.1-0.3 mags (Margot, 2005; Noll et al, 2004a,b; 
Osip et al, 2003; Grundy et al., 2007). This similarity im- 
phes that the components are composed of similar material, 
at least on the surface. It also suggests that the assumption 
of equal albedo and density usually made for binaries may 
have some basis in fact. 

Spectra are even better composition diagnostics than 
color measurements. Separate spectra of binary compo- 
nents are currently available only for the Pluto/Charon pair 
(Buie et al., 1987; Fink and DiSanti, 1988; Bide and 
Grundy, 2000) and for 2003 ELei (Barkume et at., 2006). 
Pluto and Charon have well known spectral differences that 
may be primarily related to the size threshold for retaining 
the very volatile CH4 and N2 ices found on Pluto but not on 
Charon. 2003 ELei and its larger satellite, by contrast, both 
have spectra that are dominated by water ice. 

Lightcurves are diagnostic of both compositional vari- 
ation on surfaces and of non-spherical shapes. They also 
give rotation rates. In binary systems, the rotation state of 
the components is subject to tidal evolution (see §3.5). Both 
unresolved and resolved lightcurves can be useful for ad- 
dressing these issues. 

Unresolved hghtcurves of short duration for a number of 
TNBs have been obtained, sometimes with incomplete or 
inconsistent results. Lightcurves of (47171) 1999 TC36 and 
(42355) Typhon/Echidna showed variations on the order of 
0.10-0.15 mags, but no period was determinable from the 
data (Ortiz et al, 2003). Similarly, observations of (66652) 
1999 RZ253 and 2001 QC298 showed small amphtude, but 
non-systematic variation over a 4-6 hour duration (Kern, 
2005). Romanishin et al. (2001) reported a lightcuve for the 
unresolved binary (26308) 1998 SMies, obtained from the 
1.8m Vatican Advanced Technology Telescope in 1999 and 
2000, with a moderate amplitude of 0.56 mags. The period 
was determined to be either 3.983 hours (single-peaked, 
caused by an albedo spot) or 7.966 hours (double-peaked, 
caused by nonspherical shape). The single-peaked period is 
near the break-up period (3.3 hours) for a solid ice body. 
Because the unresolved lightcurve of (26308) 1998 SMies 
did not show any color variation with time Romanishin et 
al. argued for the longer, double-peaked period as the most 
hkely. Lightcurve measurements of the same binary made 
at Lowell Observatory in 2006 found a slightly longer pe- 
riod of 8.40±0.05 hours (Spencer et al, 2006). 

Resolved ground based lightcurves of binaries are chal- 
lenging and can only be obtained under excellent condi- 
tions at a few facilities, and only for objects with the widest 
separations. Discovery observations of the binary (88611) 
2001 QT297 at Las Campanas Observatory with Magellan 
revealed brightness changes in the secondary component of 



0.3 mag in 30 minutes (Osip et al, 2003). Follow-up ob- 
servations revealed the secondary to have a single-peaked 
period of 5-7 hours while the magnitude of the primary 
remained constant. Additional resolved color lightcurve 
measurements found the two surfaces to share similar col- 
ors throughout the rotation indicating homogeneous, simi- 
lar surfaces (Kern, 2005). Similar observations showed that 
both components of 2003 QY90 to be variable. The pri- 
mary and secondary components were observed to change 
by 0.34±0.12 and 0.9±0.36 mags, respectively, over six 
hours of observation (Kern and Elliot, 2006). The large 
amplitude of the secondary component sometimes results in 
the secondary being brighter than the primary. Both com- 
ponents of the wide binary 2005 EO304 are variable with 
variations on the order of 0.3 mags over a period of 4 hours 
(Kern, 2005). 

Space based observations from HST resolve the compo- 
nents of binaries and can constrain the variabihty of com- 
ponents in these systems. The best studied system, by far, 
is the Pluto/Charon binary where detailed lightcurve mea- 
surements have been made (Buie et al, 1997). TNBs that 
have had their orbits measured by HST have multiple-epoch 
photometric measurements, although frequently the tempo- 
ral sampling is poor. (58534) Logos/Zoe shows variability 
in the primary of at least ~ 0.8 mag, making it challeng- 
ing at times to distinguish the primary from the secondary 
(Noll et al, 2004). However, with only a few widely spaced 
samples, this remains, for the moment, only an intriguing 
suggestion of a lightcurve. Three objects, (47171) 1999 
TC36, 2001 QC298 and (65489) Ceto/Phorcys have virtually 
no variation in flux, implying they may be relatively spher- 
ical, homogeneous, and/or pole-on (although we note the 
contradictory ground-based observations of (47171) 1999 
TCae). Once again, the sampling density is far to small to 
allow anything more than informed speculation. 

3.7 Orbit Plane and Mutual Events 

The tremendous scientific benefit that can derive from 
mutual occultations or echpses between a poorly-resolved 
object and its satellite was abundantly illustrated by the se- 
ries of mutual events between Pluto and Charon during the 
1980s (Binzel and Hubbard, 1997). As discussed before, 
these events enabled measurement of the sizes and albedos 
of both objects, of their distinct surface compositions, and 
even of albedo patterns on their surfaces. 

For observable mutual events to happen either the ob- 
server or the Sun (or both) must be temporarily ahgned with 
a TNB's orbit plane. An "occultation-type" event occurs 
when one component of the TNB passes in front of, and 
fully or partially occults, the other component from the ob- 
server's point of view. An "eclipse-type" event takes place 
when the TNB components are aligned with the Sun and 
the shadow of one falls on the other. Because the Sun and 
Earth have nearly equal fines of sight to TNBs, most mu- 
tual events observable from the Earth are combinations of 
occultation-type and eclipse-type events. 
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The larger the objects are compared with their separa- 
tion, the farther the orbit plane can deviate from either of 
these two types of alignments and still produce an observ- 
able mutual event. The criteria for both types of events can 
be expressed as Ri + R2 > s sin(0) where Ri and R2 
are previously defined and, ,s is their separation during a 
conjunction (equal to the semimajor axis, for the case of a 
circular orbit), and (p is the angle between the observer or 
the Sun and the plane of the binary orbit. During any con- 
junction when either criterion is satisfied, a mutual event 
can be observed. The period during which the orbit plane is 
aligned closely enough to the Sun's or to the Earth's line of 
sight to satisfy the criteria for events can be thought of as a 
mutual event season. 

Each orbit of a transneptunian secondary brings a su- 
perior conjunction (when the primary is closer to the ob- 
server) and an inferior conjunction (when the secondary is 
closer), so shorter orbital periods lead to more frequent con- 
junctions and associated opportunities for mutual events. 
The most recent mutual event season of Pluto and Charon 
lasted from 1985 through 1990, and, since their mutual or- 
bit period is only 6.4 days, there were hundreds of observ- 
able events during that season. For more widely separated, 
smaller pairs, with longer orbital periods, the mutual event 
seasons may be shorter and conjunctions may be less fre- 
quent, leading to far fewer observable events, or even none 
at all. For example, (26308) 1998 SM165 has a reasonably 
well-determined orbit with a period of 130 days, a semima- 
jor axis of 1 1,170 km , and an eccentricity of 0.47 (Margot, 
2004). From Spitzer thermal observations, the diameters of 
the primary and secondary bodies are estimated to be 294 
and 96 km, respectively (Spencer et al, 2006). Ignoring un- 
certainties in the current orbital elements, the mutual event 
season will last from 2020 through 2026, with 12 mutual 
events being observable at solar elongations of 90 degrees 
or more. Of these 12 events, 2 are purely occultation-type 
events and one is purely an eclipse-type event. The rest in- 
volve combinations of both occultation and eclipse. 

4. BINARY FORMATION 

When the Pluto/Charon binary was the only example of 
a true binary in the solar system (true binary = two objects 
orbiting a barycenter located outside either body) it could 
be discounted as just another of the peculiarities associated 
with this yet-to-be-dwarf planet. The discovery of numer- 
ous similar systems among TNOs, however, has changed 
this calculus. Any successful model for producing TNBs 
cannot rely on low -probability events, but must instead em- 
ploy processes that were commonplace in the portion of the 
preplanetary nebula where these objects were formed. For- 
mation models must also account for the observed prop- 
erties of TNBs including the prevalence of similar-size bi- 
naries, the range of orbital eccentricities, and the steeply 
increasing fraction of binaries at small angular separations. 
Survival of binaries, once they are formed, is another impor- 



tant factor that must be considered when, for instance, com- 
paring the fraction of binaries found in different dynamical 
populations or their distribution as a fraction of Hill radius. 

Several possible modes for the formation of solar sys- 
tem binaries have been discussed in the literature includ- 
ing fission, dynamical capture, and collision (c.f. reviews by 
Richardson and Walsh, 2006; Dobrovolskis et al., 1997). 
For TNOs, capture and/or collision models have been the 
most thoroughly investigated. Fission is unlikely to be im- 
portant for objects as large as the currently known TNBs. 
Other possible mechanisms for producing binaries, e.g. 
volatile-driven splitting, as is observed in comets, have not 
been explored. Interestingly, both capture and collisional 
formation models share the requirement that the number 
of objects in the primordial Kuiper Belt (at least the small 
ones) be at least a couple of orders of magnitude higher than 
currently found in transneptunian space. It follows that all 
of the TNBs observed today are primordial. 

4.1 Capture 

Capture models rely, in one form or another, on three- 
body interactions to remove angular momentum and pro- 
duce a bound pair. As we show in detail below in §4.5, they 
are also very sensitive to the assumed velocity distribution 
of planetesimals. Goldreich et al. (2002, hereafter G02) de- 
scribed two variations of the three-body model, L^, involv- 
ing three discrete bodies, and L^s, where the third body is 
replaced by a dynamical drag coefficient corresponding to a 
"sea" of weakly interacting smaller bodies. In G02 's analy- 
sis, the L^s channel was more efficient at forming binaries 
by roughly an order of magnitude. 

Astakhov et al. (2005) extended the capture model by 
exploring how a weakly and temporarily bound (a ^ Rn) 
pair of big bodies hardens when a third small body ("in- 
truder") passes within the Hill radius (O, of the pair Their 
calculations assume the existence of transitory binaries that 
can complete up to ^10 mutual orbits before the third body 
approaches. They find that a binary hardens most effec- 
tively when the intruder mass is a few percent that of a big 
body (this result likely depends on their assumed approach 
velocities of up to Suh where vn = JIk^h is the Hill ve- 
locity, with — 27r/(200yr) denoting the local Kepler 
frequency). 

The two capture formation channels and L^s require 
that binaries form — and formation times increase with de- 
creasing separation a — before v > v-a- Given the pre- 
ponderance of classical binaries {Stephens and Noll, 2006; 
Fig.), we can speculate that the primordial classical belt 
may have enjoyed dynamically cold (sub-Hill or marginally 
Hill) conditions for a longer duration than the primordial 
scattered population. 

Dynamical capture is the only viable formation scenario 
for many TNBs with high angular momentum (§4.4) and 
is a possible formation scenario for most, if not all, known 
TNBs. Given the apparent importance of capture models. 
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we explore them in detail in §4.5 with a specific focus on 
the case of binaries with similar mass components. 

4.2 Collision 

ColUsional models were proposed for the Pluto/Charon 
binary early on based on the angular momentum of the 
system (McKinnon, 1984, 1989) and one such model has 
recently been shown to be numerically feasible {Canup, 
2005). However, as we note below, angular momentum ar- 
guments alone are not sufficient to prove an impact origin. 
Nix and Hydra, the outermost satellites of Pluto, plausibly 
resulted from the same impact that generated Charon {Stern 
et ai, 2006). Ward and Canup (2006) propose that colli- 
sional debris within the exterior 4: 1 and 6: 1 resonances of 
Charon — resonances stabilized by Charon's initially large 
eccentricity — accumulated to form Nix and Hydra. Accord- 
ing to their scenario, as Charon's orbit tidally expanded, the 
small satellites would have been forced outwards to their 
current locations in resonant lockstep with Charon. Tidal 
circularization of Charon's orbit would have eventually re- 
leased Nix and Hydra from resonance. In this scenario, the 
nearly circular orbits of Nix and Hydra result from their 
coalescence from a dissipative, nearly circular disk; their 
eccentricities are not altered by resonant migration because 
the resonances involved are of the co-rotation type. 

The small satellites of Pluto-sized TNOs 2003 ELgi and 
(136199) Eris, characterized by satellite-to-primary mass 
ratios of ~1%, might also have formed by impacts (Stern, 
2002). Some collisional simulations (Durda et al., 2004; 
Canup, 2005) reproduce such low mass ratios. Tidal expan- 
sion of satellite orbits can explain, to within factors of a few, 
the current semimajor axes of the companions of 2003 ELgi 
(Brown et al. 2005a). An unresolved issue is the origin of 
the small, but significant, orbital eccentricity, 0.05 ± 0.003, 
for the outermost satellite of 2003 ELgi. Tides should have 
reduced the eccentricity to values much smaller. Also the 
mutual orbital inclination of the satellites of 2003 ELgi, 
which might be as large as 39° (Brown et al. 2006a), has yet 
to be explained. Formation by coUisionless capture along 
the lines of G02, though not ruled out, remains poorly ex- 
plored for unequal mass components (Brown et al. 2005a). 

To occur with reasonable frequency, collisions between 
Pluto-sized (R ^ 1000 km) TNOs must be gravitationally 
focused. Transneptunian space may have been populated by 
a few dozen such objects (Kenyon and Luu 1999). If their 
relative velocities were less than the Hill velocity, uh, then 
the collision timescale would be ^--^6 Myr ( fTSl ). Otherwise 
the collision time exceeds ^500 Myr Like coUisionless 
capture (see §4.1), binary formation by giant impacts must 
have taken place while the disk was dynamically cold. 

4.3 Hybrids 

Hybrid collision/capture models are possible as well; 
two variants on this theme have been proposed. 



Weidenschilling (2002) considered a model in which 
a third big body collides with one member of the scat- 
tering pair Since physical collisions have smaller cross- 
sections than gravitational interactions, this mechanism re- 
quires ^10^ more big (R ~ 100 km) bodies than are cur- 
rently observed to operate at the same rate as L^. It also 
predicts an unobserved prevalence of widely separated bi- 
naries. 

Funato et al. (2004) proposed that observed binaries 
form by the exchange reaction Ls + L L"^ + s wherein a 
small body of mass Msm, originally orbiting a big body of 
mass A/big, is ejected by a second big body. In the majority 
of ejections, the small body's energy increases by its orbital 
binding energy ^Msmvlsd'^-' leaving the big bodies bound 
with separation a ^ (Mbig/Msin)i?. As formulated, this 
model predicts a prevalence of very high eccentricity bina- 
ries that is not observed. The rate-limiting step for the ex- 
change reaction model is the formation of the pre-existing 
(Ls) binary, which requires two big bodies to collide and 
fragment ( fTsT l. 

4.4 Angular Momentum 

It is possible, in some cases, to distinguish between cap- 
ture and collision based on the angular momentum of the 
binary. Early theoretical arguments in favor of a collisional 
origin for Pluto/Charon were based on the fact that the total 
angular momentum of the system exceeds breakup for a sin- 
gle, reconstituted object (e.g. McKinnon, 1989). However, 
it can also be shown that to have formed from a fragmen- 
tary collision, binary components cannot have too much 
angular momentum. It is conventional to express angu- 
lar momentum as J / J' , where the combined orbital and 
spin angular momentum of the binary J is normalized by 
J' = \J GM^Q^Rcs, where G is the gravitational constant, 
Aftot is the total system mass, and Rcs is the radius of 
an equivalent spherical object containing the total system 
mass. Canup (2005) found that binary systems produced 
by single collisions have J/ J' < 0.8 (for an order-of- 
magnitude derivation of this result, see Chiang et al., 2006). 
For instance, the Earth/Moon system has J/ J' ~ 0.1 and 
the Pluto/Charon system has J/ J' ~ 0.4. 

In Table 2, we list J/ J' for TNBs. We use actual mea- 
surements where they exist, otherwise the calculation as- 
sumes that both binary components have spin periods of 8 
hours and densities of 1 g cm^''. For about half of the TNBs 
documented in Table 2, J/ J' exceeds unity, so much angu- 
lar momentum that formation via two-body collisions can 
be ruled out. 

4.5 Detailed Capture Models 

Capture models are of particular importance for TNBs 
as the only class of models that can explain the existence of 
high angular momentum systems. We review and expand, 
in detail, on capture models in this section. 
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TNOs can become bound ("fuse") by purely gravita- 
tional means while they are still dynamically cold. Fol- 
lowing G02, we consider how "big" TNOs, having sizes 
^big ~ 100 km, fuse when immersed in a primordial sea 
of "small" bodies, each of size i?sm- We assume that the 
small bodies contain the bulk of the disk mass: the surface 
density of small bodies a ^ ctmmsn, where ctmmsn ~ 0.2 
g cm^^ is the surface density of solids in the minimum- 
mass solar nebula at a heliocentric distance of 30 AU. We 
assume that the surface density of big bodies was the same 
then as it is now: S ^ O.OIctmmsn- This last condition 
agrees with the output of numerical simulations of coagula- 
tion by Kenyan and Luu (1999). The velocity dispersion 
of small bodies is u > uh. For convenience we define 
a = Rhig/Rn — 1-5 x 10^^ and note that the surface 
escape velocity from a big body Wosc — v^a^^^^. The ve- 
locity dispersion of big bodies is v < u. 

Small bodies have their random velocities u amplified by 
gravitational stirring by big bodies and damped by inelastic 
collisions with other small bodies. Balancing stirring with 
damping sets u {G02): 



fRs 



1/4 



v-1/2 



Rs 



20 m 



1/4 



(10) 



Big bodies have their random velocities v amplified by 
gravitational stirring by other big bodies and damped by 
dynamical friction with small bodies. When v > wh, this 
balance yields (Goldreich et al. 2004, hereafter 004): 



u \<y ) 3 
Combining ( fTOl ) with (fTTT ). we have 



(11) 



1/2 



1/4 



20 m 



1/4 



(12) 

valid for _Rsm > 20 m. If i?sin < 20 m, then v < uh, 
neither dTTJ nor ([T2]i holds, but ([TOll still does. 

By allowing for the possibility that v > wh, we depart 
from G02. When v > wh, inclinations and eccentricities of 
the big bodies' heliocentric orbits can be comparable (G04). 
If, prior to fusing, big bodies have an isotropic velocity 
dispersion, then the resultant mutual binary orbits will be 
randomly inclined, in agreement with observation. Other- 
wise, if w < uh, big bodies collapse into a vertically thin 
disk (G04) and mutual orbit normals, unless subsequently 
torqued, will be parallel, contrary to observation. Invoking 
V > vu comes at a cost: Efficiencies for fusing decrease 
with increasing v. To quantify this cost, we define a nor- 
malized velocity parameter F as 



1 if V < vn 



V /uh if u > wh 



(13) 



and derive how the rates of fusing depend on F. How F 
increased to its current value of ^10"^ — i.e., how the Kuiper 
belt was dynamically heated — remains contested (Chiang 
et al., 2006; Levison et al., 2006). 

Both the and L^s scenarios described by G02 be- 
gin when one big body {L) enters a second big body's {L) 
sphere of influence. This sphere has radius R\ ~ Ru_/F^. 
Per big body, the entry rate is 



^^^.a-^F-\ 



pRhig 



(14) 



If no other body participates in the interaction, the two big 
bodies pass through their spheres of influence in a time 
^cnc ~ Ri/v ^ ^K^^ (assuming they do not collide). 
The two bodies fuse if they transfer enough energy to other 
participants during the encounter. In L^, transfer is to a 
third big body: L + L + L^L'^+L. To just bind 
the original pair, the third body must come within Ri of 
the pair The probability for this to happen in time tone is 
P]^3 Nitcnc- If the third body succeeds in approaching 
this close, the probability that two bodies fuse is estimated 
to be on the order of unity (this probability has yet to be 
precisely computed). Therefore the timescale for a given 
big body to fuse to another by is 



fuso,L3 ~ 



pRhii 



2 4 



2i^" Myr . 



(15) 

The extreme sensitivity to F in equation JTSl l impUes that 
no binaries form by once v exceeds vb_- 

In L^s, energy transfer is to smaU bodies by dynami- 



cal friction: L + L + s° 



In time tone, the 



pair of big bodies undergoing an encounter lose a fraction 
{ailK/ pRhig){vcsc/u)'^tcnc of their energy, under the as- 
sumption Vcsc > u > wh {G04). This fraction is on the 
order of the probability ^^2^ that they fuse, whence 



^fusc. L^s 



1 



Rsi 



1 



20 m 



pi?big 



Myr , 



2 p 2 
Rhig ^i- 



-F' 



(16) 



where we have used (fTOl i. The steep dependence on F im- 
plies that formation by L^s remains viable only for u/uh 
less than a few. Together equations ( fTSl l and ( fT6] l imply 
that explaining random binary inclinations by appealing to 
V >VYi spawns a fine-tuning problem: Why should w « wh 
during binary formation? 

Having formed with semimajor axis a ^ Ri ^ 
(7000/F^)i?big, the mutual orbit shrinks by further en- 
ergy transfer If is the more efficient formation process, 
passing big bodies predominantly harden the binary; if L^s 
is more efficient, dynamical friction dominates hardening. 
The probability P per orbit that a shrinks from r^Ri to 
^Ri/2 is of order either P^a or Pl^s- We equate the for- 
mation rate of binaries, A^aii/^fusc, with the shrinkage rate. 



17 



flKPNbin\x~Ri, to conclude that the steady-state fraction 
of TNOs that are binaries with separation Ri is 



/bin (a ^ -Ri) 



bin I 



all 



pRhig 



OAF- 



(17) 

As a decreases below Ri, shrinkage slows. Therefore /bin 
increases with decreasing a. Scaling relations can be de- 
rived by arguments similar to those above. If L^s domi- 
nates, /bin a" for a < Ruivu/u)"^ and /bin oc for 
a > Ruivn/u)'^ (G02). If dominates, /bin oc a-^/^ 
when V < vh and /bin oc when v > vu- For reference, 
resolved TNBs typically have a ~ lOOiJbig- The candidate 
close binary reported by Sheppard and Jewitt (2004), and 
any similar objects, may be the hardened end-products of 

and L^s (though a coUisional origin cannot be ruled out 
since for these binaries J/ J' might be less than unity). 

These values for /bin(«) characterize the primordial 
disk. Physical colhsions with small bodies over the age 
of the solar system, even in today's rarefied environment, 
can disrupt a binary. For this reason. Petit and Mousis 
(2004) find that the widest and least massive binaries, hav- 
ing a > 400i?big, may originally have been ten times more 
numerous than they are today. 

Using the formahsm we have developed for capture, it 
is also possible to develop expressions for the collision 
timescale. As with capture, this rate of collisions is sensi- 
tive to the relative velocity of potential impactors, F. This 
results from the dependence of the colUsional rate on grav- 
itational focusing. For the two different velocity regimes, 
the timescales are given by 



^fuse, exchange '"^ * 



^a3/2 _ 0.6Myr if v < vu 



(18) 



5. THE FUTURE OF BINARIES 



Binaries offer unique advantages for the study of the 
Kuiper Belt and are hkely to be among the most intensively 
studied TNOs. Statistical studies will continue to refine the 
range of properties of binaries, their separations, their rel- 
ative sizes, and their frequency as a function of dynamical 
class, size, and other physical variables. The frequency of 
multiple systems will be constrained. As binary orbits are 
measured, orbital parameters such as semimajor axis, ec- 
centricity, and orbit plane orientation will also become the 
subject of statistical investigations. System masses derived 
from orbits will fuel vigorous studies of physical proper- 
ties of objects, both singly and in ensemble. Albedo and 
density are constrained by the measurement of system mass 
alone, and can be separated with the addition of thermal in- 
frared measurements. The internal structure of TNOs can 
be inferred from densities. The extremely low densities 



of a few objects measured so far requires structural mod- 
els with a high fraction of internal void space. The study 
of lightcurves of binaries is at a very early stage and can 
be expected to shed light on the shapes, pole orientations, 
and tidal evolution of binaries. This, in turn, may yield ad- 
ditional information on the internal structure of TNOs by 
constraining the possible values of Q. 

There are probably more than 100 binary systems that 
could be discovered in the currently-known transneptunian 
population of 600+ objects with well-estabhshed heliocen- 
tric orbits. As the TNO population expands, the number of 
binaries can be expected to expand with it. Most of these 
discoveries will be made with HST or other instruments 
with equivalent capabilities. There are also yet-to-be ex- 
plored areas of interest where theory and modeling can be 
expected to make significant progress. As progress is made 
in understanding how binaries formed in the Sun's proto- 
planetary disk, these principles can be extended to other cir- 
cumstellar disks that are now found in abundance. If binary 
protoplanets are common, as seems to be the case for the 
solar system, we may even expect to find binary planets as 
we explore extrasolar planetary systems. 
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